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Abstract 



For large classes of systems of polynomial nonlinear PDEs necessary 
and sufficient conditions are given for the existence of solutions which 
are discontinuous across hyper-surfaces. These PDEs contain the 
Navier-Stokes equations, as well as those of General Relativity and 
Magneto- Hydro dynamics . 



1. Preliminary Remarks 

There has for longer been an interest in finding junction condition 
across hyper-surfaces of discontinuities for solutions of various non- 
linear equations of Physics, such as in General Relativity, Magneto- 
Hydrodynamics, and so on. Two approaches in this regard have been 
pursued in the literature. One of them is trying to keep the equations 
and introduce hard to deal with integral formulations of the respective 
conditions, while the other is introducing ad-hoc simplifying assump- 
tions which are dictated by mathematical convenience, rather than 
physical reasons. 

In view of the obvious disadvantages of both such approaches, in 
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Braunss, a direct method was suggested for a certain limited class of 
nonlinear equations, aiming to obtain the desired junction conditions 
from the formulation of weak solution conditions across the respective 
discontinuities. 

However, in Braunss, this approach faces the following problem which 
is not addressed in a rigorous way in view of the fact that the equa- 
tions dealt with are nonlinear. Namely, due to the discontinuities 
across the hyper-surfaces, the weak solutions include Heaviside func- 
tions. Thus in the respective equations one becomes involved with 
nonlinear operations, among them products with Heaviside functions 
and their derivatives, the Dirac distributions, as well as possibly with 
the derivatives of Dirac distributions. Consequently, the respective 
weak solution conditions fail to have a rigorous formulation in terms 
of the customary and essentially linear theories of Schwartz distribu- 
tions or Sobolev spaces. 

Furthermore, the class of nonlinear PDEs dealt with in Braunss is 
rather limited. 

This direct method of Braunss in formulating the junction conditions, 
a method with an obviously sound and worthwhile underlying idea, 
was approached rigorously for the first time in Rosinger [6, pp. 139- 
162], see also Rosinger [7,8], where it was applied to the large class of 
arbitrary polynomial nonlinear PDEs, which thus contain among oth- 
ers the Navier-Stokes equations, as well as those of General Relativity 
and Magneto-Hydrodynamics. 

Such a rigorous approach was made possible by the use of the differen- 
tial algebras of generalized functions first introduced and developed in 
Rosinger [1-15], see also Mallios & Rosinger [1-3], Rosinger & Walus 
[1,2], Mallios [1,2]. 

One of the consequences is the identification of a very large class of 
nonlinear systems of PDEs, called resoluble, for which such junctions 
conditions exist. 

Furthermore, an explicit characterization for the existence of such 
junction conditions is obtained for that class of nonlinear systems of 
PDEs. 

Resolution of singularities in the above context means, therefore, the 
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rigorous expression of the conditions required upon solutions in the 
mentioned algebras across hyper-surfaces of discontinuities. 

The main results in this regard are reviewed here, and the respective 
proofs can be found in Rosinger [6-8]. 

Certain improvements and clarifications in the formulation of the ear- 
lier results are presented. 

Related to the above, one should note the following less well known, 
yet nontrivial fact. Certain possible restrictive conditions imposed on 
the classes of nonlinear systems of PDEs may happen to come from 
Physics. Indeed, as shown by Rubel for ODEs, and by Buck for PDEs, 
and mentioned briefly in section 10, there are such equations which 
may seem to have in a certain sense too many solutions. 



2. The Class of Polynomial Nonlinear Systems of PDEs 

The polynomial nonlinear systems of PDEs dealt with are of the form 



(2-1) Ei 



<i<h 



a 



Cp,i( X ) 111 



<i<*.,, U a ^ (x) 



where x G f2 C ]R n , with Q nonvoid, open. The unknown functions 
are U = (E/i, ... , U a ) : Q. — > R a , while 1 < ap tiJ < a, 1 < (3 < 
b, Pp,i,j € N n , and the functions cg ; j, fp £ C°°(fi) are given. 



The system of PDEs (2.1) is called of type (ME), if and only if each 
of the associated partial differential operators 



(2.2) T p (x,D)U(x) = E {cp,i( x ) [ill- r /;/ '" ' U api 

can be written in the form 



(x 



']} 



Tp(x,D)U(x) = 

(2.3) 

= Ei< P < r/3 Lp ;P (x,D) [ U asp (x)Pf,, p {x,D) U ap p (x) ] 

where Lp jP (x,D) are linear partial differential operators with C°°- 
smooth coefficients and of order mp, p , while Pp }P (x,D) are linear 
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partial differential operators with C°°-smooth coefficients and of or- 
der at most one. 

It is easy to see that the Navier-Stokes equations, as well as the equa- 
tions of General Relativity and those of Magneto-Hydrodynamics are 
of such type (MH). 

The condition (2.3) can be written in the equivalent form where 

u «e, P (x)Pp, p (x,D) U ap p (x) 

are replaced respectively with 

(2-4) Ei<„,„'< fl U a (x) Pp, p , a , a ' (x,D) U a , (x) 

here P/3 ;P , a ,a' ( X ,D) being again linear partial differential operators 
with C°°-smooth coefficients and of order at most one. 

3. The Singularities 

We turn now to the class of singularities dealt with. As is well known, 
for every closed subset r C f2, there exists a C°°-smooth function 
7 : Q — > R, such that 

(3.1) T = { x G VI | 7 (a:) = } 

Here however, we shall impose the restriction that 

(3.2) T has zero Lebesgue measure 

4. The Type of Discontinuous Solutions 

The solutions U = (U u ... ,U a ):Sl — > R a of (MH) type systems of 
nonlinear PDEs (2.1) we are interested in are supposed to satisfy 

(4.1) U are C°° - smooth on Q\T 
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thus they can only have discontinuities across the singularity sets T. 
Consequently, we shall be looking for such solutions which are of the 
form 

(4.2) U(x) = U-(x) + [ U + (x) - U_(x) } H 7 (x), xett 

where U- and U + are classical solutions of the PDEs (2.1) on the whole 
of Q, while the Heaviside type functions : Q — > R are defined by 



(4.3) H y (x) 



if j(x) < 

1 if 7(ar) > 



5. Formulation of the Problem 

We are looking for the necessary and sufficient junction conditions 
satisfied by 

• the (MH) type nonlinear systems of PDEs (2.1), 

• the singularities Y and their representations 7 in (3.1), 

• the classical solutions U- and U + of the PDEs (2.1) on the whole 
of Q, 

so that U in (4.2) are solutions of the PDEs (2.1) on the whole of f2. 

6. Characterization of the Existence of Junction Conditions 

In suitable differential algebras of generalized functions, one obtains 
Theorem 6.1. 

The function U in (4.2) is a global solution on Q of the (MH) type 
nonlinear system of PDEs (2.1), if an only if the following junction 
conditions are satisfied for each 1 < j3 < b, in the neighbourhood of 
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the singularity set T, namely 

E p L PtP { [£ Q>Q , [(U + ) a P[3, p , a , a ,(U + ) al - (U-) a Pf l , p ,a,a>(U-) a ,]] H,) + 
+ \ Z P L(3, P { [E a , a - [(U + ) a + (U-U [{U + ) a> - (U-) a ,]] Q(3, p , a , a >H,} = 

where Q p, Pta ,a' is the first order homogeneous part of P/3, P;a ,a'- 
Remark 6.1. 

In view of the way the Heaviside type functions H 1 appear in the 
junction conditions above, it is obvious that these junction conditions 
are trivially satisfied outside of the singularity sets T, that is, on the 
open and dense subsets fl \ T. 

Therefore, these junction conditions act only in a neighbourhood of 
the singularity sets T. 

The fact that their actions are not restricted to the singularity sets V 
alone, but also involve their neigbourhood, results from the structure 
of the differential algebras of generalized functions within which these 
junction conditions hold. 

It should be noted that such an involvement of a neighbourhood of 
singularities, and the impossibility to reduce the issues involved only 
to the singularity sets is typical for various customary theories of gen- 
eralized functions, and in particular, for the Schwartz distributions or 
the Sobolev spaces. For instance, in the case of the one dimensional 
Dirac 5 distribution, this distribution vanishes on R \ {0}. However, 
its behaviour at G R can only be described if a neighbourhood of 
that point is considered, and not only that point alone. 

7. Resoluble Systems of Nonlinear PDEs 

The above results can further be extended to a yet more large class of 
nonlinear systems of PDEs given by 

Definition 7.1. 
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The polynomial nonlinear system of PDEs (2.1) is called resoluble, if 
and only if the respective nonlinear partial differential operators (2.2) 
are of the form 

Tp{x,D)U{x) = 

(7.1) 

for iGfl and 1 < j3 < b, whenever 

(7.2) U(x) = 4)(x) + x(x)u(x), xeVl 

with C°°-smooth x '■ ^ — > ^ a an d oj — > R, while pp, p G 
N ra , lp tP EN, and Tp jP (x, D) are arbitrary C°°-smooth coefficient poly- 
nomial nonlinear partial differential operators in -0 and X- 

Proposition 7.1. 

Every (MH) type nonlinear system of PDEs is resoluble. 

□ 

The main result, which extends Theorem 6.1., is the following 
Theorem 7.1. 

The function U in (4.2) is a global solution on f2 of the resoluble 
nonlinear system of PDEs (2.1), if an only if the following junction 
conditions are satisfied for each 1 < (3 < b, in the neighbourhood of 
the singularity set T, namely 

(7-3) Ei< P < r/3 [Tp, P (x,D) (U-,U + -U-)]D>*pH, = f fi 

□ 

The junction conditions (7.3) in Theorem 7.1. contain as a particular 
case those in Theorem 6.1., and obviously have a more simple, com- 
pact form. In the next section, indications are given on the way these 
junction condition (7.3) can effectively be computed. 
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8. Computation of the Junction Conditions 

The relations which are useful for the effective computation of the 
junction conditions in Theorem 7.1. are the following 

(8.1) D p H^ = (D p j)5^ P eN n , \p\ = 1 

(8.2) DP(D'<J 7 ) = ( J D p 7 )(D' +1 5 7 ), p G N n , |p| = 1, leN 

(8.3) 7 (5 7 = 

(8.4) 7j D' +1 <5 7 + (Z + 1)L>'5 7 = 0, leN 

where 5 1 is the Dirac distribution concentrated on the hyper-surface 
T. In this regard, it is well known that (8.1), (8.2) are valid when 

(8.5) grad^f(x) ^ 0, x G T 

In view of (8.1), (8.2), one easily obtains 

(8.6) ZM 7 = Eo<i<| P |-i[^P,i(^^)7]^7. PGN»,| P |>1 

where K p j(x,D) are polynomial nonlinear partial differential opera- 
tors with C°°-smooth coefficients and order at most / + 1, which can 
be obtained from the recurrent relations 

(8.7) Kp Q^XjD)^ = -D P 7, p G N", | p | = 1 

Kp+^D^ = D^Kp^D^} + [K Ptl . 1 {x,D) 1 \D*' 1 , 

(8.8) 

p,q eN n , | q | = 1 

Now in view of (8.6), the junction conditions (7.3) can further be sim- 
plified. 

9. Reduced Smoothness 
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As shown in Rosinger [6, pp. 121-162], the above results hold for the 
more general case of polynomial nonlinear systems of PDEs (2.1) when 
the coefficients are merely continuous. In such however, one 

has to make use not of a suitable single differential algebra of gener- 
alized functions, but instead, of chains of such algebras. 



10. Too Many ODEs and PDEs ? 

In Rubel it was shown that there exist nontrivial fourth order polyno- 
mial nonlinear ODEs 

(10.1) P(U(x),U'(x),U"(x),U"'(x),U""(x)) = 0, i£l 

where P is an integer coefficient polynomial in four variables, with the 
following dense solution property. 

Given any two continuous functions /, to G C°(R), with w(x) > 0, for 
i£l, then there exists a solution U G C°°(R) of (10.1), such that 

(10.2) | f(x) - U(x) | < w(x), xeR 

Two further related facts can be mentioned. The ODEs (10.1) can be 
constructed in a rather simple manner. However, they are implicit in 
the highest derivative U "". 

A corresponding version of the above dense solution property was ob- 
tained for PDEs in Buck. Namely, for every n G N, n > 2, there exist 
nontrivial polynomial nonlinear PDEs 

P(U(x), ... ,DPU(x), ...) = 0, 

(10.3) 

p G N n , \p \ < m(n), x G [0, l] n 

with P a real coefficient polynomial of degree at most d(n), such that 
for every continuous function / G C°( [0, l] n ) and e > 0, there exists 
a solution U G C°°( [0, 1 f) of (10.3), with the property 

(10.4) | f(x) - U(x) | < e, x G [0, l] n 



9 



This time, however, the proof is rather involved, as it uses Kolmogorov's 
solution to Hilbert's Thirteenth Problem. Also, the degrees d(n) of 
the polynomials P, and the orders m{n) of the PDEs (10.3) are very 
large, even for small values of n, such as n = 2, 3, 4. 
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